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We compute exactly the mean perimeter and area of the convex hull of N independent planar 
Brownian paths each of duration T, both for open and closed paths. We show that the mean 
perimeter (Ljv) = oln VT and the mean area (An) — PnT for all T. The prefactors cvn and /3n, 
computed exactly for all N, increase very slowly (logarithmically) with increasing N. This slow 
growth is a consequence of extreme value statistics and has interesting implications in ecological 
context in estimating the home range of a herd of animals with population size N. 

PACS numbers: 05.40.-a, 02.50.-r, 87.23.Cc 



Ecologists often need to estimate the home range of an 
animal or a group of animals, in particular for habitat- 
conservation planning [1]. Home range of a group of ani- 
mals simply means the two dimensional space over which 
they typically move around in search of food. There ex- 
ist various methods to estimate this home range, based 
on the monitoring of the positions of the animals over a 
certain period of time [2j] . One method consists in draw- 
ing the minimum convex polygon enclosing all monitored 
positions, called the convex hull. While this may seem 
simple minded, it remains, under certain circumstances, 
the best way to proceed [3[. The monitored positions, for 
one animal, will appear as the vertices of a path whose 
statistical properties will depend on the type of motion 
the animal is performing. In particular, during phases 
of food searching known as foraging, the monitored po- 
sitions can be described as the vertices of a random walk 
in the plane |4|, [5j . For animals whose daily motion con- 
sists mainly in foraging, quantities of interest about their 
home range, such as its perimeter and area, can be es- 
timated through the average perimeter and area of the 
convex hull of the corresponding random walk (Fig. 1 1(a)] ). 
If the recorded positions are numerous (which might re- 
sult from a very fine and/or long monitoring), the num- 
ber of steps of the random walker becomes large and to a 
good approximation the trajectory of a discrete-time pla- 
nar random walk (with finite variance of the step sizes) 
can be replaced by a continuous-time planar Brownian 
motion of a certain duration T. 

The home range of a single animal can thus be char- 
acterized by the mean perimeter and area of the convex 
hull of a planar Brownian motion of duration T start- 
ing at origin O. Both 'open' (where the endpoint of the 
path is free) and 'closed' paths (that are constrained to 
return to the origin in time T) are of interest. The latter 
corresponds, for instance, to an animal returning every 
night to its nest after spending the day foraging in the 
surroundings. For an 'open' path, the mean perimeter 
(Li) = v / 87rTandthe mean area (Ai) = ttT/2 are known 
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FIG. 1: Convex hull of (a) a 7-step random walk starting at 
O (b) 3 closed Brownian paths. 



in the mathematics literature H, 0, 
only the mean perimeter is known 



For a 'closed' path, 
: (Li) 



In any given habitat, an animal is however hardly alone 
but they live in herds with typically a large population 
size N. To study their global home range via the con- 
vex hull model, one needs to study the convex hull of TV 
planar Brownian motions. Assuming that the animals 
do not develop any substantial interaction between them 
during foraging, this leads us to the main question ad- 
dressed in this Letter: what is the mean perimeter and 
area of the convex hull of N independent planar Brown- 
ian motions (both open and closed) each of duration T? 
This question is of vital importance in ecological conser- 
vation: if the population size increases by, say ten-fold, 
by how much should one increase the home range, i.e., 
the conservation area? 

Using the standard scaling property of Brownian mo- 
tion, lengthscale ~ (timescale) 1 / 2 , it follows that the 
mean perimeter and area of the global convex hull of 
N independent Brownian paths will behave as (Ln) = 
oinVT and (An) = PnT for all T. The main challenge is 
to estimate the N- dependence of the prefactors a at and 
Pn • The central result of this Letter is to provide exact 
formulae for a at and /3n for all N, both for open and 
closed paths. 
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Most interestingly we find that the mean perimeter 
and area of the convex hull increases very slowly with N 
for large N: a at — 27rVzln7V, /3n — 27rlniV for open 
paths, and a/v(c) ~ iwjl In iV, /3n(c) 77 IniV for closed 
paths, c referring to closed paths. This leads to our main 
conclusion: the home range increases very slowly with 
increasing population size. Indeed, the In N behaviour of 
the prefactor (3n(c) m the mean area indicates that, for 
instance, a 10-fold increase in the number of animals in 
the herd will result, on average, in the addition of only 
about 3.6 T units of area to the home range of the herd 
— a good news for conservation. 

To proceed, let I denote any set of points on the plane 
with coordinates (x^yi). Let C denote the convex hull 
of /, i.e., the minimal convex polygon enclosing this set. 
Geometrically, the convex hull can be constructed very 




FIG. 2: Support function of a 7-step random walk 

simply: consider any arbitrary direction from the origin 
O specified by the angle 9 with respect to the x axis. 
Bring a straight line from infinity perpendicularly along 
direction and stop when it touches a point of the set I 
[e.g., in Fig. [21 this point has co-ordinates (xk* , Vk*)]- Let 
M{6) denote the Euclidean distance of this perpendicular 
line from the origin when it stops, measuring the maximal 
extension of the set / along the direction 0. This support 
function can be simply written as 

M{0) = max {xi cos + yi sin 0} . (1) 

Knowing M(9), one can express the perimeter and the 
area of the convex hull C via Cauchy's formulae [lQj : 
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(M 2 (0) - (M'(9)) 2 ) , (3) 



where M'{0) = dM/dO. In this Letter we show how these 
two formulae can be used to compute the mean perime- 
ter and area of the convex hull of N planar Brownian 
motions. 




FIG. 3: (a) Location r* of the maximum M{0) of zg(r) and 
(b) corresponding value of M'(0) = ho(r*) 



To illustrate the main idea let us start with a single 
open Brownian path, the generalizations to TV > 1 and 
for closed paths will follow. In this case the set / consists 
of the vertices of a Brownian path of duration T start- 
ing at the origin O. Since it is a continuous path, we can 
conveniently label the coordinates by B(r) = (x(t), y(r)) 
with < r < T. Here x(r) and y{r) are just two indepen- 
dent one dimensional Brownian paths each of duration T 
and evolve via the Langevin equations: x(r) = rj x (r) 
and y(r) = rj y (r) where rj x (r) and rj y (r) are independent 
Gaussian white noises, each with zero mean and delta 
correlated, e.g., {t] x {t)t] x {t')) = S(r — r') and the same 
for rj y . Clearly then (x 2 (r)) = r and {y 2 (r)) = r. 

Let us now consider a fixed direction 0. Then, zq(t) = 
x(r) cos Q + y{r) sin 9 and he (r) = —x(r) sin + y(r) cos 
are also two independent one dimensional Brownian mo- 
tions (each of duration T) parametrized by 6. Then M{6) 
in Eq. (pQ) is simply the maximum of the one dimensional 
Brownian motion zq(t) over the time interval r G [0,T], 
i.e., M{0) = max[^(r)]. Let r* be the time at which 
this maximum is achieved. Then, M(ff) — zq(t*) = 
x(t*) cos# + y(r*) sin 9. Taking derivative with respect 
to gives M'(0) = -£(Y*)sin<9 + cos 6> = h e (r*). 
Thus, while M{6) is the maximum value of the first 
Brownian motion zq{t\ M'(0) is the value of the sec- 
ond independent Brownian motion he (r) but at the time 
r = r* at which the first one achieves its maximum (see 
Fig. 3(a) and |3(b)[ ). Note, in particular, that for = 0, 
zq(t) = x{t) and ho(r) = y(r). Thus M(0) is the max- 
imum of x(t) in r G [0, T] while M'(0) = y(r*) is the 
value of y at the time r* when x achieves its maximum. 

Taking expectation and using isotropy of the mean over 
the choice of directions, it follows that 



(L) = 2tt (M (0)) 
(A) = ^«[M(0)] 2 > 



<[M'(0)] 2 » . 



(4) 
(5) 



For a one dimensional Brownian motion x(r) over [0,T], 
the distribution of its maximum is well known [11]. The 
cumulative distribution, Qi(m,T) = Prob[M(0) < m] = 
erf(ra/>/2T), where erf(z) = ^ f* e'^ du. The first 
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two moments can be easily computed: (M (0)) = y /f 2T/ir 
and ([M(0)] 2 ) = T. Eq. (|4]) then gives the mean perime- 
ter, (Li) = aiVT with ai = V^Stt- The calculation of the 
mean area is slightly trickier since we need ([M'(0)] 2 ). 
We first note that for any fixed time r*, E[y 2 {r*)\ = r* 
(since it is a free Brownian motion) where the expectation 
is taken over all realizations of the process y at fixed r* . 
But r* itself is a random variable, being the time at which 
the first process x achieves its maximum. Its distribution 
is also well known and given by the celebrated arcsine law 
of Levy [12]: P(r*,T) = [r*(T - r*)]- 1 / 2 /^. Thus, av- 
eraging further over r* taken from this distribution, we 
finally get ([M'(0)] 2 ) = (r*> = T/2. Eq. © then gives 
the exact mean area (A\) = f3\T with f3\ = tt/2. 

For a single closed path, the analysis is similar, ex- 
cept that x(r) and y(r) are now two independent one 
dimensional Brownian bridges of duration T, i.e., both 
start at the origin but are constrained to return to the 
origin at time T: x(0) = x(T) = and y(0) = y(T) = 0. 
The distribution of the maximum of a Brownian bridge is 
also well known and its first two moments are given 
by: (M (0)) = y^^T/S and ([M(0)] 2 ) = T/2. Thus, the 
mean perimeter of the convex hull is given from Eq. (|4|): 
(Li) = ai(c)VT with a\{c) = ^tt 3 /2. To compute the 
mean area, we note that for a Brownian bridge y(r), at 
fixed time r*, E[y 2 (r*)} = r*(T - r*)/T. Moreover, the 
distribution of r* , the time at which a bridge achieves its 
maximum is known to be uniform [11]. Thus, taking av- 
erage over r* with uniform distribution, P(r*, T) = 1/T, 
we get (M'(0) 2 ) = (?/ 2 (t*)) = T/6. Eq. © then gives 
the exact mean area of the convex hull of a planar closed 
Brownian path: (A^ = /?i(c)T with ^(c) = tt/3. To 
our knowledge, this result, even for a single closed path, 
appears to be new. 

This method can then be generalized to N independent 
planar Brownian paths, open or closed. We now have 
two sets of N Brownian paths: Xj(r) and jjj(r) (j = 
1,2, ...,7V). All paths are independent of each other. 
Since isotropy holds, we can still use Eqs. (|4]) and ([5j), 
except that M(0) now denotes the global maximum of 
a set of N independent one dimensional Brownian paths 
(or bridges for closed paths) Xj(r) (j = 1,2,..., TV), each 
of duration T, 

M(0) = max max \x\(r). ^(r), . . . , xn(t)] . (6) 

1<3<N 0<r<T L V n V n V /J V 7 

Let j* and r* denote the label of the path and the time 
at which this global maximum is achieved. Then, using 
argument similar to the N = 1 case, it is easy to see that 
M'(0) = yj m (r*), i.e., the position of the j*-th y path at 
the time when the x paths achieve their global maximum. 

To compute the first two moments of M(0), we first 
compute the distribution P/v[M(0), T] of the global max- 
imum of N independent Brownian paths (or bridges) 
Xj (r) . This is a standard extreme value calculation. Con- 
sider first N open Brownian paths. It is easier to compute 



the cumulative probability, Qjv(ra,T) = Prob[M(0) < 
m). Since the Brownian paths are independent, it fol- 
lows that Q N (m,T) [Qi(ra, T)]^, where Qi(m,T) = 
erf(m/v // 2T) for a single path mentioned before. Know- 
ing this cumulative distribution Qat(M(0), T), the first 
two moments (M(0)) and ([M(0)] 2 ) can be computed for 
all TV. Using the result for (M(0)) in Eq. (|4]) gives us 
the mean perimeter, (Ln) = ajyVT with 

ql 7v = 47V / duu e~ u2 [erf (it)]*' 1 . (7) 
Jo 

The first few values are: a\ — = 5.013.., ct2 — 

4V7T = 7.089.., a 3 = 24 tan^l/V^)/^ = 8.333.. etc. 
(see Fig. 4 for a plot of a at vs. N). For large TV, one can 
analyse the integral in Eq. by the saddle point method 
giving [13j], djy — 2i\\/2 ln N. This logarithmic depen- 
dence on N is thus a direct consequence of extreme value 
statistics [l4| and the calculation of the mean perimeter 
of the convex hull of N paths is a nice application of the 
extreme value statistics. 

To compute the mean area, we need to calculate 
([M^O)] 2 ) in Eq. ©. We proceed as in the N = 1 
case. For a fixed label j and fixed time r, the expecta- 
tion E[y 2 (r)] = r which follows from the fact that yj(r) 
is simply a Brownian motion. Thus, E[y 2 ^{r*)] = r*. 
Next, we need to average over r* which is the time at 
which the global maximum in Eq. (|6]) happens. The 
distribution P/v(r*, T) of the time of global maximum of 
N independent Brownian motions, to our knowledge, is 
not known in the probability literature. We were able to 
compute this exactly for all N [13]. Skipping details, we 
find that P N (r*,T) = ±F n (t*/T) where 

/>OC 

F N (z) = - g . / dxxe-^ [erf^)]^ 1 (8) 
V 2(1 - z ) Jo 

and cin is a normalization constant fixed by, 
Jq Fn(z) dz = 1. It is easy to check that for N = 1, it 
reproduces the arcsine law mentioned before. Averaging 
over r* drawn from this distribution, we can then 
compute ([M^O)] 2 ) = / Q T r* P N (r*, T) dr\ Substituting 
this in Eq. ([5]) gives the exact mean area for all TV, 
(A N ) = p N T with 

rOO 

f3 N = 4N^ J duu [erf^)]^" 1 (ue~ u2 - g(u)J (9) 
where q(u) = ttt= fn e / ' dt • For example, the first few 

/ 2^F JO y/t{l-t) F ' 

values are given by, f3\ = ir/2 — 1.570.., $2 = tt = 3.141.., 
p 3 = tt + 3-V^ 4.409.. etc. (Fig. 4 shows a plot of (3 N 
vs N). The large TV analysis gives |13| . (3n — 27rlniV. 
Thus for large TV, the shape of the convex hull approaches 
a circle of radius a/2 ln N which, incidentally, coincides 
with the set of distinct sites visited by N Brownian mo- 
tions fist- 



For N closed Brownian planar paths one proceeds in 
a similar way. Without repeating the analysis, we just 
mention our main results [13[. The mean perimeter and 
area are given by, (Ln) = aN(c)VT and (An) — /3n(c)T 
where, for all TV, 
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and f(k) = Q (k - I)- 3 / 2 (k tan-^Vife 
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The first few values are: ai(c) = y/7r 3 /2 = 3.937.., 
a 2 (c) = Vtt3(a/2- 1/2) = 5.090.., a 3 (c) = Vt^(S/V2 - 
3/2 + I/a/6) = 5.732.. and /?i(c) = tt/3 = 1.047.., 
/?2(c) = tt(4 + 3tt)/24 = 1.757.., /3 3 (c) = 2.250.. etc. 
(see Fig. 4 for a plot of aAr(c) and /?jv(c) vs. A 7 "). 
Large iV analysis shows that [13[, a/v(c) c± 7rV2 In TV and 
Pn(c) — § lniV, smaller respectively by a factor 1/2 and 
1/4 than the corresponding results for open paths. 

We have also computed the mean perimeter and area 
of the convex hull of N = 1, 2, 3, 4 Brownian paths (both 
open and closed) via numerical simulations. For each TV, 
we constructed TV independent normal Gaussian random 
walks of 10 4 steps each with a time step At = 10 -4 . 
For each realisation of the TV walks, we constructed the 
convex hull using the Graham scan algorithm [16] and 
computed its perimeter and area and then averaged over 
10 3 samples. We find excellent agreement with our ana- 
lytical predictions (Fig. Hj). 

The method presented here is general and can, in prin- 
ciple, be applied to compute the mean perimeter and area 
of the convex hull of any set of random points. In fact, 
when the set consists of just random points, each drawn 
independently from a given distribution, the statistics 
of the perimeter and area of their convex hull has been 
studied extensively in various different contexts [17[ . Our 
method, using extreme value statistics, can easily recover 
these results [13j, but can go further e.g., when the points 
are correlated as in the case of Brownian paths. 

This work leads to several interesting open questions. 
It would be interesting to extend the results presented 
here for normal diffusion to the case where the trajec- 
tories of animals such as birds, deers or spider mon- 
keys undergo anomalous diffusion, e.g., Levy flights 
etc. 



18 



Il9l . |20[ . Another interesting question concerns 
the effect of interactions or collective behavior of the an- 
imals on the statistics of their convex hulls. For animals 
like birds that move in 3-dimensional space, it would be 
interesting to study the statistics of the convex polytope 
of their trajectories, such as the mean surface area and 
the volume of such a polytope 0. Finally, the distri- 
bution of the time of maximum r*, a crucial ingredient 
in our method, has recently been computed exactly for 
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FIG. 4: Setting T = 1, the analytical results for average 
perimeter q>n [Eq. J7J)] and area [3n [Eq. ©] of N open 
(Op.) Brownian paths, and similarly the average perimeter 
a N (c) [Eq. (HUJ)] and area (3 N (c) [Eq. JTTJ)] of N closed (CI.) 
Brownian paths, plotted against N. The symbols denote re- 
sults from numerical simulations (up to N = 4) . 



a variety of constrained one dimensional Brownian mo- 
tions [2l[ . These results may be useful to study the statis- 
tics of convex hulls of constrained planar Brownian paths 
(see e.g. [22|). 

We thank D. Dhar and H. Larralde for useful discus- 
sions. 
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